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Abstract. We justify an applicability of the adiabatic perturbation theory for three well known 
systems with impacts: a ball between two slowly moving walls, a slowly irregular waveguide, and 
an adiabatic piston. 

The adiabatic perturbation theory (see, for example, [T|, p. 200 - 211) is used to describe 
the dynamics in smooth Hamiltonian systems containing variables of two types: fast and 
slow. The accuracy estimates of this theory are known. In a number of cases, one can for¬ 
mally apply the procedure of this theory to discontinuous Hamiltonian systems, in particular 
to systems with impacts. However, the validity of such formal approach does not follow from 
the theory for smooth systems. A distinctive feature of considered problems is fast (instant) 
variation of “slow” variables at an impact. In the present work we obtain accuracy estimates 
of the adiabatic perturbation theory for systems with impacts. We consider three model 
problems that belong to the three main classes of systems, where the adiabatic perturbation 
theory is used: 1) the motion of a ball between two slowly moving walls corresponds to the 
Hamiltonian systems depending on slowly varying parameters; 2) the problem of a slowly 
irregular waveguide corresponds to the Hamiltonian systems slowly depending on some of 
variables ; 3) the problem of an adiabatic piston corresponds to the Hamiltonian systems 
with slow and fast motions. We discuss these problems in order to demonstrate the effec¬ 
tiveness of the approach based on the adiabatic perturbation theory: the dynamics in these 
problems is well known, and it was described earlier by other methods (in particular, in [2] 
for problem 3)). 

The first approximation of the procedure of the adiabatic perturbation theory leads to 
a conclusion that the system has an adiabatic invariant (an approximate integral). Such 
a conclusion is often used for systems with impacts, but its validity is derived from direct 
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calculations (see ^j, Sec. 52). Higher approximations of the procedure of the adiabatic 
perturbation theory for systems with impacts are formally considered in [I]. The perturbation 
theory for non-smooth Hamiltonian systems is discussed in works yet only in the 

case of continuous phase variables. In systems with impacts some of phase variables are 
discontinuous at a point of reflection. 

Dealing with systems with impacts, one usually considers Poincare map in order to obtain 
results analogous to those obtained for smooth systems by means of higher approximations 
of the adiabatic perturbation theory. If this map is smooth, one can either use the procedure 
of the adiabatic perturbation theory for smooth maps (see, for example, £Zj), or use an 
artificial approach, namely, to describe the map as a Poincare map for an auxiliary smooth 
Hamiltonian system and apply the procedure of the adiabatic perturbation theory to this 
system 0. 

An advantage of the approach described in the present paper is that one can deal with 
Hamiltonians of original systems, and consider systems with impacts similarly to smooth 
systems. It is no need to construct and deal with Poincare maps (that can turn out to 
be non-smooth as, for example, in the problem of an adiabatic piston, sec.4). Under this 
approach calculations become simpler and more concise. 

1 Fermi-Ulam model 

The problem of oscillations of a particle between two slowly moving walls is called the Fermi- 
Ulam problem (or model) (FigIH)^J. Let the left wall be fixed (this assumption is made to 
shorten calculations) and the right wall slowly change its position. The distance between 
the walls is d(et) > const > 0, where t is the time variable, £ > 0 is a small parameter. 
Introduce slow time variable r , t — et . Let function d(-) G C°° . We consider the problem 
on a slow time interval that either does not depend on £ or grows as £ decrease. In the 
second case we shall suppose that function d and its derivatives are bounded on the real 
axis. Let the mass of the particle be equal to 1, then the velocity of the particle is equal to 
its momentum. One can describe the dynamics of the particle in the following way: between 
the walls the particle has a constant velocity v ; the velocity of the particle change sign after 
an impact with the left wall, and transforms as v\ = 2d — v after an impact with the right 
moving wall. The Hamiltonian of the particle is well defined between walls and describes 
free motion of the particle with a constant momentum v . 

One can interpret the problem as the motion of a particle in a potential that equals 0 
between the walls and +oo in the other space. The right wall is moving slowly. Therefore, 
it is useful to consider the problem when the right wall is fixed (d = const). The phase 
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portrait of this system is shown in FigJ21 On the portrait one can define “action-angle” 
variables (/, 0) in a standard way. “Action” is the area bounded by a phase trajectory 
divided by 2tt : I — -E ■ \v\ ■ 2d — ^v / 2 E, where E is the Hamiltonian of the particle. 

2 r2 

Hence, E = ^4- . The angle (phase) 0 is the uniformly changing angle variable on the 
trajectory; 0 = 2i r|r, where t is the time interval corresponding to the motion of the 
particle from the initial point to a given point and T is the period of motion. 

Suppose the phase is zero at the left wall; then 


*1’ 


0 G (0, 7r), 


tt(2-§), 0g(tt, 2tt), 


( 1 . 1 ) 


where x is the distance between the particle and the left wall. 

If the phase is defined by cud, then generating function W = W{x, /, d) of the form 


W = 


7r 


lx 
d ’ 


( 1 . 2 ) 


0 g (0, tt), 

tt(2 - f ), </> e (t 2tt) 

transforms variables (v, x ) to the “action-angle” variables. 

Use generating function W to make a change of variables in the real system with right 
wall moving slowly. The motion of the particle between the walls is described by the 


Hamiltonian system with Hamiltonian H = E + ^ : 


TT 2 / 2 Id 

H= u?- ~d m ’ 


(1.3) 


where 


m 


0, 0 e (0, tt), 

0 — 27r, 0 G (7T, 27r). 


(1.4) 


Hamiltonian 03 completely describes the motion of a particle between the walls. 

At an impact with the fixed wall, the value of the “action” variable is preserved. The 
variation of the “action” variable after an impact with the moving wall is 


U -1- = 


2 dd 


7r 


(1.5) 


where /_ and I + are the values of / prior to and after the impact correspondingly. 

Let us describe an impact by means of Hamiltonian 03- An interaction with the wall 
is instantaneous. Hence, it is likely that the actual value of the variation of the “action” 
variable after an impact with the moving wall can be calculated in the following way: fix 
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the moment of time at the impact, consider Hamiltonian ra as a Hamiltonian of a one 
d.o.f. Hamiltonian system, and obtain the required value of the “action” variable using the 
energy conservation law. 

Lemma 1 This approach properly describes jump of the “action” variable at an impact with 
the moving wall. 


Proof. 

Suppose the value of time variable is fixed. Set the value of energy prior to an impact 
equal to the value of energy after an impact: 


n 2 P nl_d 7r 2 / 2 | irl + d 

~2d 2 cT = + d 

Hence, we obtain 


- ip 

2d 




from which immediately follows 03. 

As applied to the considering problem, the basic statement of the adiabatic perturbation 
theory can be formulated as follows. 


Theorem 1 After canonical change of variables (/, (jf) 


(/, 0) with generating function 


11 — If) + £S(I , 0, T, £), S — Si + £^2 + . . . + £ r 2 SV_1, 

where r is any fixed natural number, Hamiltonian in takes the form 


( 1 . 6 ) 


'H = r, e) + e r H r (i, 0(J, 0, r, e), r, e), 


(1.7) 


He, r = E + eH i + ... + £ r " 1 K r _ i . 

Here Si — Si(I,<f>,T) , Hi = Tt-fl , r) . Functions Si, Hi are C°° in I , r; functions Si 
are continuous with respect to 0. 


Proof. 

We use the standard procedure of the adiabatic perturbation theory (see, for example, 
P). Make a canonical change of variables with generating function (H~6l) in the system with 
Hamiltonian 03. The variables transform according to the following expressions 

(1- 
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/ = / + £ 


dS_ 

00 ’ 


— 0 + £- 


dS_ 

~di ‘ 







New Hamiltonian 7 i has the form 


Tt = 77(7, 0, t, e) + e 2 


dS(I,</>,T,e) 

dr 


(1.9) 


Define functions S t so that the new Hamiltonian has the form O- Substitute expres¬ 
sions 03 in EH) and equal terms of the same order in e . Thus, we find the set of equations 
from which S'*. Ti t can be found. For example, in the first order in e we get 


Ui (7,r) = 


dE 3S 1 


T'» 


9 / 50 

A prime denotes derivative with respect to r. After averaging over 0 one obtains 
Tii — 0. Therefore, Si can be found as a quadrature, which turns out to be continuous 
with respect to 0. One can choose Si so that its average is 0. Similarly, one can derive 
functions Si that are continuous with respect to 0 and C°° with respect to I , r. The 
required properties of functions 7 % , H r follow. 


Corollary 1.1 The value of variable I is preserved along a solutioii with an accuracy 0(e r ) 
on time interval 0(e ~ k ) for any prefixed natural number k . The value of variable 0 is de¬ 
fined by an integrable Hamiltonian system with Hamiltonian 7Ts.r with an accuracy 0(e r ~ k ) 
on the same time interval. 


Proof. 

Consider the change of variables (7,0) i—► (7,0) generated by function W = 70 + 
eS(1, 0, r, e), S = Si + eS 2 + ... + £ r+k ~ 2 S r „i . 

According to Theorem d the Hamiltonian function in the new variables has the form: 


H 


We, r+kp, t, e) + £ r+t H r+ t(I,<l>, t, e). 


( 1 , 10 ) 


The expressions for the change of variables give 1 — 1 = 0(e r ) . Hence, below we consider 
variation of variables 7 and 0. Variation of value 7 between the impacts with the walls 
is 0(e r ). Consider variation of value 7 at an impact with the moving wall. Lemma H) 
describes the way to calculate the value of variable 7 after an impact with the wall. In 
order to determine variation of variable 7 one substitutes 7 = 7(7, 0, r, e) in Hamiltonian 
77 in (11.31) and notes that 77 = Ti — £ 2 (P . Thus, one obtains the following expression 


7f S)r+fc (7 + , r, e) + e r+k H r+k (I + , tt + 0, r, e) - 


2 dS(I + , 7T,T,e) 


dr 


— 77s,r+fc(7-, T , e) + 


t> 





~\~S 


r +k h 


r-\-k 


(/_, 7T — 0, r, e) — £ 


2 dS(I_,7T,T ,£) 


<9r 


Here 7_ and 7 + are the values of 7 prior to and after the impact respectively. We note 
that the function S is continuous with respect to 0 ■ Hence, we find 


d E(L,r ) 
dl 


+ o( £ 2 ) I (/+-/-) = 0(0. 


Therefore 7 + — 7_ = 0(V +fc ). 

The particle collides with the walls not more than 0(e k ) times. Therefore the value of 
7 varies as 0(e r ). Hence the part of the statement concerning variable 7 is proved. 

Consider now variation of 0 . On any time interval between the impacts with the moving 

A A ‘3'iy A 

wall value of 0 deviates from the solution of equation 0 = V r , 7 = const by 0(e r ). At 
a collision with the moving wall value of variable 0 jumps. Value of the jump can be found 
from the formulas for change of variables CD and depends on variation of 7. Functions 
S m are continuous with respect to 0. Hence, the jump in the phase at an impact with the 
wall is 0(e r+l ). This completes the proof of the statement. 


Corollary 1.2 Value of variable I is preserved with an accuracy 0(e) on a time interval 
0(s ~ k ) for any prefixed natural number k . 

Corollary 1.3 The formulas for change of variables 7,0 i—> 7,0 and approximate expres- 

a /v Q'kL 

sions I = const, 0 = ^ describe the behaviour of variables I,<f> with an accuracy 0(£ r ) 

for I and 0(s r ~ k ) for 0 on the time interval 0(s~ k ) . 

Remark 1 ZemmaQ theorem^ and corollaries \1.1\ - fO can be used in case of motion 
in potential U(x,t ) G C°° between the walls, and the particle in adiabatic approximation 
collides with the walls at a nonzero velocity. 

Remark 2 In I7W systems with impacts are investigated in the following way. Instead of 
system with impacts, a smooth system with a repulsive potential strong in a thin layer (the 
thickness of the layer S << 1, the gradient of the potential ~ 1/5 ) is considered. Properties 
of the system with impacts are derived from that of the smooth system proceeding to limit as 
5 —> 0 . Following this approach one could define such a potential for Fermi-Ulam problem 
(and for problems considered below, sec. 2,3), use the perturbation theory for the smooth 
system and obtain corresponding accuracy estimates for the system with impacts. However, 
one should consider uniformity of accuracy estimates in 5 . The method defined above allows 
to avoid this problem. 
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2 A slowly irregular planar waveguide 


In this section we consider slowly irregular planar waveguides. The problem is to find the 
trace for a ray of light propagating in a slowly irregular waveguide with reflective walls. 
A slow irregularity means that the width of the waveguide is changing slowly along the 
waveguide’s length EH- 

Assume, for simplicity, that one wall of the waveguide coincides with x axis. The 
position of another wall in the plane Oxy is defined by formula y = d(X ), where X — ex. 
e is a small parameter (FigCH). Let d(-) G C°° . Let us consider the problem on a time 
interval 0(e ~ k ), where k is any prefixed natural number. Assume that function d, its 
derivatives of any order and function 1/d are uniformly bounded on the real axis. 

One can describe the propagation of rays of light in a medium using the Hamiltonian 
system with Hamiltonian 


H = pl+p 2 y -n 2 (x,y), (2.1) 

where one should consider only zero energy level H = 0 EH- Here p XtV are the variables 
canonically conjugated to the coordinates x,y and n(x,y ) is the refraction index. In our 
case inside the waveguide n — 1. The Hamiltonian system is not well defined at the reflective 
walls and one can use only conservation laws there. In the unperturbed system X = const 
and the projection of a phase trajectory on p y , y plane looks similar to the phase trajectory 
of a particle in FigEl Like in the previous section, one can make a canonical change of 
variables (p y , y) to “action-angle” variables (/, 0): / = Py ^ /1A:> ; if phase 0 is zero at the 
bottom wall, then 



0 G (0, 7r), 


0 G (7T, 27t). 


( 2 . 2 ) 


Using standard definition of the phase JH one obtains the generating function for the 
change of variables (p y , y ) to the “action-angle” variables in the following form 


W(I,y,X) 


7T^, 0G(O, 7T), 

7r ( 2_; ?)’ 0e(vr, 2vr). 


(2.3) 


Make a canonical change of variables ( p x , x . p y , y) i—> ( p x , x , / , 0 ) with generating 
function S = p x x + W (/, y, ex) in the exact (perturbed) system. The Hamiltonian in the 
new variables is 


H 


7T 2 / 2 

~/P~ 


Px 


Id> f(A s 
£ —- fw 


1, 


(2.4) 
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where is defined in (II.dj) . the prime denotes a derivative with respect to A". 

At a reflection at the bottom wall value of the “action” is preserved. At a reflection at 
the upper wall it transforms according to the following expression 


h 


I - 


‘led! 

1 + e 2 d 


£ d'I + 



(2.5) 


Remark 3 Value p x is preserved at a reflection. Formula (HH) defines the law of trans¬ 
formation for the angle of the ray of light arriving to the bottom wall. This angle transforms 
after an impact with the upper wall as a\ = a — 2 arctg{ed'), thus implying \2.A) . 

Lemma 2 Consider Hamiltonian \2.f | ) and fix variables (p x , x ) . Value of the Hamiltonian 
is preserved when value of </> passes through i r . This law of conservation defines the actual 
value of “action’ variable I after reflection at the upper wall of the rvaveguide. 

Proof. 

Set equal values of Hamiltonian EH) before and after reflection of the ray at the upper 
wall: 


7T 2 / 2 

d 2 


Id 1 , 1 

Px ~ E ~d' K ~~ 1 


vrV , (, , hd' , , 

— + { p * +£ —* ] - L 


After transformations obtain 


7T 2 / 2 2 j/2 \ „ „ Id' -K 2 1 i 2 ,/2\ 0 „ h d' 

— —{1+e d ) - 2ep x —7T = — — (1 + e d ) + 2e P x ~r 77 - 


As value of the “action” variable is positive, we find 


h = I 


‘led! p x d 


1 + e 2 d! 7T 


Substituting the expression for momentum p x = p x — £^j-7r , we finally get formula 

r t 2ed ' ( W t d P*\ 

h ^ , -iT^ {ed, + — ) ’ 

coinciding with EH). The lemma is proved. 

Let us formulate the basic statement of the adiabatic perturbation theory for the con¬ 
sidered problem. 

Theorem 2 After canonical change of variables ( p x ,x,I,cj )) i—► (p x ,x,I,(p) with generat¬ 
ing function 


W = I(/> + p x x + eS(I,f>,p x ,ex,e), 


S = S 1 + eS 2 + ... + e r ~ 2 S^ 


( 2 . 6 ) 
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where r is any prefixed natural number, Hamiltonian takes the following form 


H = H^ r {j,p x ,ex,e) + e r H r (I, (f,p x , ex, e), 


(2.7) 


77s,r — —p -b Px — 1 + £hi\ + . . . + £ r 1 Hr-l- 

Here functions Si, TU are C°° with respect to I,p x ,ex; Si are continuous functions of 

0 - 


Proof. 

Make a canonical change of variables ( I,<f,p x ,x ) h-> (/,0,p x ,7) with generating func¬ 
tions (ESI). The variables are transformed as: 


0 = 0 + £ !?> 

Px = Px + e 2 g, 

££ = ££ + £ 2 -|^. 

OPx 


( 2 . 8 ) 


Following the standard procedure of the adiabatic perturbation theory, define functions 
S m such that the Hamiltonian in the new variables does not depend on 0 up to terms of 
order 0(e m ) and S m are continuous and periodic with respect to 0. The procedure is as 
follows: substitute formulas for the change of the “action” and longitudinal variables ( p x , x ) 
from EH) to Hamiltonian EH and choose the generating function so that the Hamiltonian 
in the new variables has form EH- 

Thus, one can recurrently define functions S m . For example, an equation for Si looks 
as follows: 


as, p x d(x)d' 

One can choose functions S m in such a way that they have zero average with respect to 
0. Finally, the Hamiltonian takes the form we need. Functions Si, Hi obviously have all 
the properties mentioned in the theorem. 

Consider the solutions of Hamiltonian equations with Hamiltonians EH and 
7 ~tz,r(J,P,£x,e) corresponding to initial conditions 7(0) = J(0), 0(0) = 0(0), p x (0) = 
p(0), 7(0) = x(0), where 0 is the phase variable conjugated to variable J. 
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Corollary 2.1 Value of variable I is preserved with an accuracy 0(e r ) on a time interval 
£~ k for any prefixed natural number k. The projection of the ray’s trajectory on p x , sx 
plane on such an interval lies in the 0{e r ) -neighbourhood of curve Hs !r (J,p,£x,£) = 0. 
Behaviour of variables £x,p x ,(f> is described by the solution of the system with Hamiltonian 
with accuracies 0(£ r ~ k+1 ), 0(£ r ~ k+1 ) , and 0 (£ r ~ k ) respectively (under some natural 
additional conditions given below). 

Proof. 

Similarly to the previous section, consider an auxiliary change of variables with generating 
function 


W = 1(f)+ p x X + £S 1 (I,(f>,p x ,£x) + . . . + £ r+k 1 S r+k _ 1 (I,(f),p x ,£x). 

According to Theorem [21 Hamiltonian (12.411 takes the form: 

H = H^ r+k (I,p x ,£X,£) + £ r+k H r+k (I,(j),p x ,£X,£), 

where 0 and x are considered as functions of new variables (I,(f>,p x ,x ). Formulas for the 
change of variables give I—I = 0 (£ r ), 0—0 = 0 (£ r ), x—x = 0 (£ r ) and p x —p x = 0(e r+1 ). 
Therefore, it is sufficient to consider behaviour of variables I,<j),p x ,x. 

Consider an approximate conservation of variable I . The ray’s Hamiltonian in the new 
variables is 


H = Hz,r+k(I,Px,ex,e) + E r+k H r+k {I,(f>,p x ,EX,£). (2.9) 

On each time interval between the reflections of the ray at the upper wall variable / vary 
according to Hamiltonian equation 


j _ c r+k 9H r+k 

<90 

After summing up over all such intervals, it gives an accuracy of conservation 0(£ r ). 

It follows from EH that values of variables (p x , x) jump at a reflection of the ray at 
the upper wall. However, the generating function for the change of variables is continuous 
with respect to 0. Therefore, formulas for the change of variables give A p x = 0(e 2 AI) , 
A ex = 0(e 2 AI) , where A/, A p x , A ex are jumps of variables /, p x . ex after reflection 
of the ray at the upper wall. Using Lemma E one estimates value of the jump of variable 
7 as 0(£ r+k ) . The ray reflects at the upper wall not more than 0 (£~ k ) times on the time 
interval £~ k . Therefore, total variation of variable I is 0 (£ r ) . Thus, the part of the 
theorem concerning variable I is proved. 


10 



Projection of the rays trajectory on p x , ex plane is in 0{e r ) -neighbourhood of the curve 
Tis t r(J,P,£x,e) = 0 because along the trajectory 7ds ir (J, p x , ex, e) = 0(e r ). 

Values of variables p, x, 0 change continuously, yet values of variables p x ,x,4> jump 
at a reflection of the ray at the upper wall. It was noted above that values of jumps 
of variables p x ,x, 0 are determined by the value of jump of I at the reflection of the 
ray at the upper wall as A p = 0(e 2 AI) = 0{e r+k+2 ) Ax = O(eAI) = 0{e r+k+l ) and 
A0 = 0[eAI) = 0(e r+k+l ). Thus, total variations of values of variables (p x , x,0) due to 
reflections are 0(e r+2 ) for variable p x and 0(e r+1 ) for variables x n 0. Total variations 
of variables p x , x , 0 due to reflections are the same. 

Consider the zero approximation of Hamiltonian 7Yv,r : 


H 


o 

S ,r 


dg p 


7T 2 / 2 

U2?-T + P 

a z {£x) 


The Hamiltonian system with Hamiltonian \'H% r describes the motion of a particle in 
potential U = 2 $lex) ■ Assume (see F igEJ) that for the ray under consideration value of 
function F satisfies one of the following conditions: 1) it is larger than value of any local 
maximum of 2 U or 2) the ray can reflect at a potential hump only once or 3) the ray moves 
between two potential humps (resonator). 

Consider case 1). The motion of the ray has a given direction, variables x and x vary 
monotonically. One can consider variables p x and p as functions of variables ex and ex 
respectively. At x = x values of variables p x and p differ by 0(e r ). It takes different 
time for variables x and x to reach a given value x* . This difference can be estimated as 
follows: 


/ dx = °( £r ~ k )• 

*(o) 

Hence, for a given moment of time t one obtains: ex(t) —ex(t) = 0(e r ~ k+1 ), p x {t) —pit) = 
0{e r ~ k+l ). 

The deviation of phases 0 — ip a t x = cc = x* is : 



u{I,p x ,ex,e) u(J,p,ex,e) 


dx + 0(e r ~ k ) = 0(e r ~ k ) 


r—k\ 


X 


X 


( 2 . 10 ) 


where we introduced u = ■ Therefore, for a given moment of time t one gets: 0(f) — 

0(f) = 0{e r - k ). 

In case 2) the ray can change its direction of propagation only once. Here one can 
consider two different regions and use either variable x or p as monotonically changing 
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variable in the corresponding region. In case 3) the waveguide is configured as a resonator. 
One can introduce an “angle” variable in the phase portrait of the system with Hamiltonian 
T-C^, r (I,p, x) as a monotonically changing variable. All arguments here are the same as in 
case 1). 

Corollary 2.2 Value of variable I is preserved with an accuracy 0(e) on a time interval 
e~ r for any prefixed natural number r . 

Corollary 2.3 Behaviour of variables I,(j),p x ,x can be found from the change of variables 
I,(j),p x ,x i—>• I,4>,p x ,x together with approximate expressions defining motion in the system 
with Hamiltonian 7 '~l^ r (I,p x ,£x) . An accuracy of such a description is the same as the 
accuracy for variables I,(f),p x ,x according to corollary 11 

Remark 4 The Hamiltonian system with Hamiltonian 

H=^+Pl-1 ( 2 . 11 ) 

considered on the energy level H = 0 describes the ray’s trajectory with an accuracy 0 (e) 
on time intervals of order e _1 . 

Remark 5 Suppose that there is nonhomogeneous medium between the walls such that a 
refraction index n(sx, y) G C°° and trace of the ray in adiabatic approximation is inclined 
to the walls at reflection at a nonzero angle; then Lemma [U Theorem^ and its corollaries 
GU-GL3 are still valid. 


3 Dynamics of a massive piston in a gas of light particles 


The problem of an adiabatic piston is an important model in statistical mechanics. It 
is considered in the context of attempts to derive thermodynamics laws from the laws of 
mechanics (see, for example, 0.D2I). The corresponding system consists of a container with 
a massive cylindrical piston and a gas of identical light particles that move independently 
elastically colliding with the walls of the container and with the piston. 

Let the mass of a particle be equal to 1. Length of the container subtracting thickness 
of the piston L and the number of gas particles are of the order 1. Mass of the piston 
M is large in comparison with mass of the gas. Suppose the piston is at rest at the initial 
moment of time; then energy of the system does not depend on mass of the piston. Hence, 
one can estimate energy of the piston as 0 ( 1 ) and its typical velocity as 0 (;^=) . Therefore, 
it is useful to introduce a small parameter e = —= . Let us discuss the problem on a time 
interval of order e _1 . 
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Below, to shorten the calculations, we consider the case when there is only one particle 
on each side of the piston (Fig0. One can easily generalize following arguments for the 
case of any prefixed number of particles. At the end of the section we give the result for the 
general case. Without loss of generality, assume that the particles’ velocities are parallel to 
the axis of the container; thus, the motion is one-dimentional. Let indexes l, r correspond to 
variables of left and right particles respectively. The variables without indexes correspond 
to the piston. Thus, full energy of the system is of the form: 


E = e 


,P 2 


pf 

+ + 


pl_ 

2 ’ 


(3.1) 


where p^ r are momenta of particles, P is momentum of the piston. Denote as x^ r distances 
between particles and the left wall of the container, as X distance between the piston and 
the left wall of the container. The distance for the right particle is calculated subtracting 
the thickness of the piston. 

Suppose the piston is fixed. For each particle make a canonical change of variables from 
(p,x) to “action-angle” variables (/, fa) the same way as in Fermi-Ulam problem (Sec. 1). 
Let phases fa )T be zero at the walls of the container; then 


vrf, 0 G(O, tt), 
vr(2-f), 0G(vr, 2vr), 


<j) r 


</>e(o, vr), 

*■(2 0€(tt, 2tt). 


(3.2) 


Conjugated “action” variables are i) = , I r = X ' 1 . If phases are defined by 

jia, generating functions for the change of variables (p, x ) to “action-angle” have the 
following forms: 


Si(h,x b X) 


S r (I r , x r , X) 


K Ij f -, (h e (0, tt), 

7 rA( 2 -f), fa e (vr, 2 tt), 

n Ir L^x r) ’ & G (°> 7r )> 

7r/ r (2 - j^Zx), fa e (7T, 27 t). 


Make the change of variables with generating function W = PX + Si(Ii,xi,X) + 
S r (I r , x r , X) in the exact system, when the piston is moving. The Hamiltonian in the 
new variables has the form: 


n = £ 2 - P 


^fifa) + L ^ x f{<t>r 


7 PI 2 

+ TTTTTr + 


2 t2 


7 PI, 


2X 2 2(L — X) 2 ’ 


(3.3) 


where function f(fa) defined by (II .41) . 
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Let us find jump of the “action” variable at an impact. It is natural to fix values of 
all variables except the variables of the colliding particle, obtain the Hamiltonian for the 
particle and consider the conservation law for the value of the Hamiltonian. Let us formulate 
the corresponding lemma for the left particle (the lemma for the right particle is formulated 
analogously). Let the prime denote value of the variable after an impact. 


Lemma 3 This approach defines actual value of the “action’ variable after an impact: 


i; 


2PX 2 PX 
7 tM + M + 1 TT M 


(3.4) 


Remark 6 Velocities v and V of the particles of mass 1 and M transform after an 
elastic collision as follows: 


I V' = V + J fc{v-V), 

\ v’ = 2V-v + j^(v-V). 

“Action” of the particle and momentum of the piston are Ii 
which follows \3.f\)- 


^ , P = MV , from 


Remark 7 If two particles collide with the piston at the same moment, the further dynamics 
is not defined. Measure of initial conditions corresponding to such collisions equals 0 . 
Therefore, we do not consider these initial condition. 


Proof. 

An impact is the the passage of cfa through the value it . Equating values of the energy 
prior to- and after an impact, we find 


7T 2 /f M + 1 _ 7 Tl t ( p I r f (0r ) \ = TT 2 // 2 M + 1 p , l r f (<t>r) 

2X 2 M MX 1 L-X’ 2X 2 M MX 1 L-X 


(3.5) 


After calculations one obtains value of the “action” variable: 


I'i = Ii 


2 r 


M + 1 7T 


L-X 


(3.6) 


Let P denote value of momentum of the piston prior to an impact; then 


J'-J 2 X (P l 7Tl1 ) 


(3.7) 
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Finally 




and the lemma is proved. 


2PX 2 PX 

TiM + M + 1 SrM ~~ 


Theorem 3 Values of variables Ii. r are preserved with an accuracy 0(e) on a time interval 


Proof. 

Introduce the normalized momentum of the piston by expression P = eP and consider 
Hamiltonian (13.31) with an accuracy 0(e 2 ): 


H = 


P 


n 2 I? 


2 ' 2X 2 + 2(L-X) 2 +£ ^L- 


X 


f(4> r ) + 0(e 2 ). 


(3.8) 


Make a canonical change of variables ( Ii- r . 4>i,r , £ _1 P, X) i—► (Ii, r , <fi, r , £ X P , 
X ), eliminating dependence of the Hamiltonian on the phase with an accuracy e 2 . The 
generating function for this change of variables has the form: 


W 


-PX + hfa + I r <t> r + eS,(4 (t>u p, X) + eS r (I r , 0 r , P, X). 


(3.9) 


Consider the expressions connecting old and new variables: 


( 






T _ T x ®Sl,r 

71, r — 71, r T e Q J, , 

7 _ : 8Sl,r 

Vl,r 0i,r d - ^ Qp r i 


P 

X 


p c-2 9S(_ , 2dSr_ 

1 " r£ si TC ax’ 

y I r 2dSi , 2dSr_ 

A + £ dP +£ w- 


(3.10) 


Substitute these expressions in Hamiltonian (13.81) and choose functions S) ir so that the 
Hamiltonian in the new variables does not depend on the phases in the first order in e. 
Therefore, functions Si, r satisfy equations 


(Xn = PXfjfa) dSr = P(L - X)f(<j>r) 

<90/ 7T 2 ’ d(p r 7T 2 


It follows from EH) that functions Si, r are defined up to an arbitrary function of 
variables Ii, r ,P,X. Let this function be equal to 0 and call variables Ii, r the “improved 
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actions”. It is easy to see from (Em EIU), that the new phase coincides with the old one 
( <pi t r = 4>i,r ), and the Hamiltonian has a discontinuity at = n . 

Lemma 131 and formulas (18.1 011 . (18.1 II) allow to find variation of the “improved action” 
after an impact as 0(e 2 ) . Variation of the “improved actions” between collisions is also 
0(e 2 ) . One can estimate the number of collisions on a time interval e _1 as 0(e~ 1 ) . Hence, 
on a time interval £ _1 “improved actions” vary by 0(e) . The “action” variables are related 
to “improved actions” by formulas for the change of variables Ii. r = Ii jT + ■ Therefore, 

“action” variables vary by 0(e) on the same time interval and the theorem is proved. 


Corollary 3.1 Let 


2 P 2 

H = E 2 — + 


7T 


2 1? 


+ 


7r 


2 j2 


2X 2 2 (L-xy 


Ii r = const, 


(3.12) 


where I^ r are initial values of the “action’' variables. The Hamiltonian system with Hamil¬ 
tonian EH describe the behaviour of variables eP,X with an accuracy O(e) on a time 
interval E^ 1 . 


One can similarly consider the system with any prefixed number of particles. Value of 
the “action” variable of each particle is preserved with an accuracy 0(e) on a time interval 
£ -1 . The Hamiltonian system with Hamiltonian 


H = e 2 


P 2 


+ 


7T 


-T 2 + 


7T 


(3.13) 


2 M ' 2X 2 EJ 2(L — X) 2 E ’” 
describes behaviour of variables eP,X with an accuracy 0(e) on the same time interval. 
Here /1 r and If, ; are the sums of squared initial values of “action” variables for the left 
particles and for the right particles respectively. (This result was first derived in |2j by a 
different method.) 

Remark 8 In the considered approximation the piston oscillates in potential U = + 

2 (L-x)^ kr ( F wW- 
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Figure 1: A ball between slowly moving walls. Figure 2: A phase portrait for the motion of the 

particle. 




Figure 3: A planar slowly irregular waveguide. 
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Figure 5: The system of a piston and two par¬ 
ticles. 


Figure 4: The plot of function 2 U and levels of 
function F. 



Figure 6: The effective potential energy of the 
piston. 
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Abstract. We justify an applicability of the adiabatic perturbation theory for three well known 
systems with impacts: a ball between two slowly moving walls, a slowly irregular waveguide, and 
an adiabatic piston. 

The adiabatic perturbation theory (see, for example, jTj, p. 200 - 211) is used to describe 
the dynamics in smooth Hamiltonian systems containing variables of two types: fast and 
slow. The accuracy estimates of this theory are known. In a number of cases, one can formally 
apply the procedure of this theory to discontinuous Hamiltonian systems, in particular to 
systems with impacts. However, the validity of such formal approach does not follow from 
the theory for smooth systems. A distinctive feature of considered problems is fast (instant) 
variation of “slow” variables at an impact. In the present work we obtain accuracy estimates 
of the adiabatic perturbation theory for systems with impacts. We consider three model 
problems that belong to the three main classes of systems, where the adiabatic perturbation 
theory is used: 1) the motion of a ball between two slowly moving walls corresponds to the 
Hamiltonian systems depending on slowly varying parameters; 2) the problem of a slowly 
irregular waveguide corresponds to the Hamiltonian systems slowly depending on some of 
variables ; 3) the problem of an adiabatic piston corresponds to the Hamiltonian systems with 
slow and fast motions. We discuss these problems in order to demonstrate the effectiveness 
of the approach based on the adiabatic perturbation theory: the dynamics in these problems 
is well known, and it was described earlier by other methods (in particular, in |2J for problem 
3)). 

The first approximation of the procedure of the adiabatic perturbation theory leads to 
a conclusion that the system has an adiabatic invariant (an approximate integral). Such 
a conclusion is often used for systems with impacts, but its validity is derived from direct 
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calculations (see [2], Sec. 52). Higher approximations of the procedure of the adiabatic per¬ 
turbation theory for systems with impacts are formally considered in [3]. The perturbation 
theory for non-smooth Hamiltonian systems is discussed in works M, yet only in the 
case of continuous phase variables. In systems with impacts some of phase variables are 
discontinuous at a point of reflection. 

Dealing with systems with impacts, one usually considers Poincare map in order to obtain 
results analogous to those obtained for smooth systems by means of higher approximations 
of the adiabatic perturbation theory. If this map is smooth, one can either use the procedure 
of the adiabatic perturbation theory for smooth maps (see, for example, 0) , or use an 
artificial approach, namely, to describe the map as a Poincare map for an auxiliary smooth 
Hamiltonian system and apply the procedure of the adiabatic perturbation theory to this 
system |B]. 

An advantage of the approach described in the present paper is that one can deal with 
Hamiltonians of original systems, and consider systems with impacts similarly to smooth 
systems. It is no need to construct and deal with Poincare maps (that can turn out to 
be non-smooth as, for example, in the problem of an adiabatic piston, sec.4). Under this 
approach calculations become simpler and more concise. 

1 Fermi-Ulam model 

The problem of oscillations of a particle between two slowly moving walls is called the Fermi- 
Ulam problem (or model) (FiglU) [HJ - Let the left wall be fixed (this assumption is made to 
shorten calculations) and the right wall slowly change its position. The distance between 
the walls is d(et) > const > 0, where t is the time variable, e > 0 is a small parameter. 
Introduce slow time variable r , r — et. Let function d(-) G C°° . We consider the problem 
on a slow time interval that either does not depend on e or grows as e decrease. In the 
second case we shall suppose that function d and its derivatives are bounded on the real 
axis. Let the mass of the particle be equal to 1, then the velocity of the particle is equal to 
its momentum. One can describe the dynamics of the particle in the following way: between 
the walls the particle has a constant velocity v ; the velocity of the particle change sign after 
an impact with the left wall, and transforms as v\ = 2d — v after an impact with the right 
moving wall. The Hamiltonian of the particle is well defined between walls and describes 
free motion of the particle with a constant momentum v . 

One can interpret the problem as the motion of a particle in a potential that equals 0 
between the walls and +oo in the other space. The right wall is moving slowly. Therefore, 
it is useful to consider the problem when the right wall is fixed (d = const). The phase 
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portrait of this system is shown in FigJ3 On the portrait one can define “action-angle” 
variables (/, 0) in a standard way. “Action” is the area bounded by a phase trajectory 
divided by 2ti : I = -^ ■ \v \ ■ 2d = ^y/2E, where E is the Hamiltonian of the particle. Hence, 
E = . The angle (phase) 0 is the uniformly changing angle variable on the trajectory; 

0 = 27r^ , where t is the time interval corresponding to the motion of the particle from the 
initial point to a given point and T is the period of motion. 

Suppose the phase is zero at the left wall; then 


*1’ 


0 G (0, 7r), 


tt(2-§), 0g(tt, 2tt), 


( 1 . 1 ) 


where x is the distance between the particle and the left wall. 

If the phase is defined by cud, then generating function W = W{x, /, d) of the form 


W = 


7r 


lx 
d ’ 


( 1 . 2 ) 


0 g (0,7r), 

*■(2 - f ), </> e (t 2tt) 

transforms variables (v, x) to the “action-angle” variables. 

Use generating function W to make a change of variables in the real system with right 
wall moving slowly. The motion of the particle between the walls is described by the Hamil¬ 


tonian system with Hamiltonian H = E + '4jj- : 


7T 2 / 2 Id 

H= u?- ~d m ’ 


(1.3) 


where 


m 


0, 0 e (0, 7r), 

0 — 27r, 0 G (7 r, 27 t). 


(1.4) 


Hamiltonian 03 completely describes the motion of a particle between the walls. 

At an impact with the fixed wall, the value of the “action” variable is preserved. The 
variation of the “action” variable after an impact with the moving wall is 


U -1- = 


2 dd 


7r 


(1.5) 


where /_ and I + are the values of / prior to and after the impact correspondingly. 

Let us describe an impact by means of Hamiltonian D3. An interaction with the wall 
is instantaneous. Hence, it is likely that the actual value of the variation of the “action” 
variable after an impact with the moving wall can be calculated in the following way: fix the 
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moment of time at the impact, consider Hamiltonian JO as a Hamiltonian of a one d.o.f. 
Hamiltonian system, and obtain the required value of the “action” variable using the energy 
conservation law. 


Lemma 1 This approach properly describes jump of the “action” variable at an impact with 
the moving wall. 


Proof. 

Suppose the value of time variable is fixed. Set the value of energy prior to an impact 
equal to the value of energy after an impact: 


n 2 P nl_d 7r 2 / 2 | irl + d 

~2d 2 cT = + d 

Hence, we obtain 


- ip 

2d 




from which immediately follows JO- 

As applied to the considering problem, the basic statement of the adiabatic perturbation 
theory can be formulated as follows. 


Theorem 1 After canonical change of variables (/, (jf) 


(/, 0) with generating function 


11 — If) + £S(I , 0, T, e), S — Si + £^2 + . . . + £ r 2 SV_1, 

where r is any fixed natural number, Hamiltonian in takes the form 


( 1 . 6 ) 


'H = r, e) + e r H r (i, 0(J, 0, r, e), r, e), 


(1.7) 


n E ,r = E + eH i + ... + £ r " 1 K r _i. 

Here Si — Si(I ,4 >,t) , Hi — TCi(I , r) . Functions Si, Hi are C°° in I , r; functions Si 
are continuous with respect to 0. 


Proof. 

We use the standard procedure of the adiabatic perturbation theory (see, for example, 
d]). Make a canonical change of variables with generating function (11.6(1 in the system with 
Hamiltonian CO' The variables transform according to the following expressions 

r o dS t . dS 

/ — / + £ — , 0 — 0 + £— j. ( 1 . 

<90 <9/ 
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New Hamiltonian 7 ~L has the form 


H = H(I, cp, t, e) + e 2 


dS(I,</>,T,e) 

dr 


(1.9) 


Define functions S t so that the new Hamiltonian has the form D3- Substitute expres¬ 
sions ca in JH and equal terms of the same order in £ . Thus, we find the set of equations 
from which S'*. H t can be found. For example, in the first order in £ we get 


Ui (J,r) = 


dE dS\ 


T'» 


dl d<p 

A prime denotes derivative with respect to r . After averaging over (p one obtains Hi = 
0 . Therefore, Si can be found as a quadrature, which turns out to be continuous with respect 
to (p . One can choose Si so that its average is 0 . Similarly, one can derive functions S { that 
are continuous with respect to (p and C°° with respect to I , r . The required properties of 
functions Hi, H r follow. 


Corollary 1.1 The value of variable I is preserved along a solution with an accuracy 0(e r ) 
on time interval 0(e~ k ) for any prefixed natural number k . The value of variable <p is de¬ 
fined by an integrable Hamiltonian system with Hamiltonian Hy, )T with an accuracy 0(e r ~ k ) 
on the same time interval. 


Proof. 

Consider the change of variables (/,</>) i—► generated by function W = Icp + 

eS{I, (p, r, e) , S = Si + eS 2 + ... + £ r+k ~ 2 S r „i . 

According to Theorem Cl the Hamiltonian function in the new variables has the form: 


H 


We, T+ k 0 , t, e) + £ r+t H r+t (I,<l>, t, e). 


( 1 , 10 ) 


The expressions for the change of variables give 1 — 1 — 0(e r ). Hence, below we consider 
variation of variables / and <p . Variation of value / between the impacts with the walls is 
0(e r ). Consider variation of value / at an impact with the moving wall. Lemma Q] describes 
the way to calculate the value of variable / after an impact with the wall. In order to 
determine variation of variable / one substitutes / = 1(1, (p, r, e) in Hamiltonian H in 
G3 and notes that H = H — e 2 ^ . Thus, one obtains the following expression 


Hz >r +k(I+, T e) + e r+k H r+k (I + , tt + 0, r, e) - £ : 


2 dS(I + ,lT,T ,£) 


dr 


— HY,, r +k(I-,T , £) + 


t> 





~\~S 


r+k H 


r-\-k 


(/_, 7T — 0, r, s) — £ 


2 dS(H,Tt,T,e) 


<9r 


Here 7_ and 7 + are the values of 7 prior to and after the impact respectively. We note 
that the function S is continuous with respect to 0 ■ Hence, we find 


d E(L,r ) 
dl 


+ o( £ 2 ) I (/+-/-) = 0(0. 


Therefore I + — 7_ = 0(e r+k ). 

The particle collides with the walls not more than 0(e k ) times. Therefore the value of 
7 varies as 0(e r ). Hence the part of the statement concerning variable 7 is proved. 

Consider now variation of 0 . On any time interval between the impacts with the moving 
wall value of 0 deviates from the solution of equation 0 = V r , I = const by 0{e r ). At 
a collision with the moving wall value of variable 0 jumps. Value of the jump can be found 
from the formulas for change of variables ca and depends on variation of 7. Functions 
S m are continuous with respect to 0 . Hence, the jump in the phase at an impact with the 
wall is 0(e r+l ). This completes the proof of the statement. 


Corollary 1.2 Value of variable I is preserved with an accuracy 0(e) on a time interval 
0(e ~ k ) for any prefixed natural number k . 

Corollary 1.3 The formulas for change of variables 7,0 i—> 7,0 and approximate expres¬ 
sions 7 = const, 0 = -jp- describe the behaviour of variables I,(j) with an accuracy 0(e r ) 
for I and 0(e r ~ k ) for 0 on the time interval 0(e~ k ) . 

Remark 1 Lemma[J\ theorem and corollaries \1 .1\ - (731 can be used in case of motion 
in potential U(x,t ) G C°° between the walls, and the particle in adiabatic approximation 
collides with the walls at a nonzero velocity. 

Remark 2 In fiTfy systems with impacts are investigated in the following way. Instead of 
system with impacts, a smooth system with a repulsive potential strong in a thin layer (the 
thickness of the layer S << 1, the gradient of the potential ~ 1 /5 ) is considered. Properties 
of the system with impacts are derived from that of the smooth system proceeding to limit as 
5 —> 0 . Following this approach one could define such a potential for Fermi-Ulain problem 
(and for problems considered below, sec. 2,3), use the perturbation theory for the smooth 
system and obtain corresponding accuracy estimates for the system with impacts. However, 
one should consider uniformity of accuracy estimates in 5 . The method defined above allows 
to avoid this problem. 
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2 A slowly irregular planar waveguide 


In this section we consider slowly irregular planar waveguides. The problem is to find the 
trace for a ray of light propagating in a slowly irregular waveguide with reflective walls. 
A slow irregularity means that the width of the waveguide is changing slowly along the 
waveguide’s length EH- 

Assume, for simplicity, that one wall of the waveguide coincides with x axis. The position 
of another wall in the plane Oxy is defined by formula y = d(X ), where X — ex. e is 
a small parameter (Fig El)- Let d(-) G C°° . Let us consider the problem on a time interval 
0(e ~ k ), where k is any prefixed natural number. Assume that function d, its derivatives 
of any order and function 1/d are uniformly bounded on the real axis. 

One can describe the propagation of rays of light in a medium using the Hamiltonian 
system with Hamiltonian 


H = p 2 x +p 2 y -n 2 (x,y), (2.1) 

where one should consider only zero energy level H = 0 ra. Here p X)V are the variables 
canonically conjugated to the coordinates x, y and n(x, y ) is the refraction index. In our 
case inside the waveguide n — 1. The Hamiltonian system is not well defined at the reflective 
walls and one can use only conservation laws there. In the unperturbed system X = const 
and the projection of a phase trajectory on p y , y plane looks similar to the phase trajectory 
of a particle in F igEl Like in the previous section, one can make a canonical change of 
variables (p y , y) to “action-angle” variables (/, 0): / = ^ Vy ^ d _ {X) ; if phase 0 is zero at the 
bottom wall, then 



0 g (0, tt), 


0 G (7T, 27t). 


( 2 . 2 ) 


Using standard definition of the phase EH one obtains the generating function for the 
change of variables (p y , y) to the “action-angle” variables in the following form 


W(I,y,X) 


7T^, 0G(O, 7T), 

7r ( 2_; ?)’ 0e(7T, 2vr). 


(2.3) 


Make a canonical change of variables ( p x , x , p y , y) > ( p x . x , / , 0) with generating 
function S = p x x + W (/, y, ex) in the exact (perturbed) system. The Hamiltonian in the 
new variables is 


H 


TT 2 / 2 

~/P~ 


Px 


£ —- fw 


1 , 


(2.4) 
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where is defined in (II.dj) . the prime denotes a derivative with respect to A". 

At a reflection at the bottom wall value of the “action” is preserved. At a reflection at 
the upper wall it transforms according to the following expression 


h 


I - 


‘led! 

1 + e 2 d 


£ d'I + 



(2.5) 


Remark 3 Value p x is preserved at a reflection. Formula (O defines the law of trans¬ 
formation for the angle of the ray of light arriving to the bottom wall. This angle transforms 
after an impact with the upper wall as a.\ = a — 2 arctg{ed'), thus implying \P,.h\) . 

Lemma 2 Consider Hamiltonian | ~2J^ and fix variables (p x , x ) . Value of the Hamiltonian 
is preserved when value of </> passes through n . This law of conservation defines the actual 
value of “action’ variable I after reflection at the upper wall of the waveguide. 

Proof. 

Set equal values of Hamiltonian EH) before and after reflection of the ray at the upper 
wall: 


7T 2 / 2 

d 2 


Id 1 , 1 

Px ~ E ~d' K ~~ 1 


vrV , (, , hd' , , 

— + { p * +£ —* ] - L 


After transformations obtain 


7T 2 / 2 2 j/2 \ „ „ Id' -K 2 1 i 2 ,/2\ 0 „ h d' 

——{I Fed )-2ep x —Ti =-p-{l +e 2 d ) + 2£p x —j-ir. 


As value of the “action” variable is positive, we find 


h = I 


2£d! p x d 


1 + £ 2 d! 7T 


Substituting the expression for momentum p x = p x — £ L ^-n , we finally get formula 

T T 2ed ' ( W T d P*\ 

h ^ , -iT^ {ed, + — ) ’ 

coinciding with EH). The lemma is proved. 

Let us formulate the basic statement of the adiabatic perturbation theory for the con¬ 
sidered problem. 

Theorem 2 After canonical change of variables (p x ,x,I,(j)) i—► (p x ,x,I,(p) with generat¬ 
ing function 


W = I(j> + p x X + £S(I,4>,p x ,£X,£), 


S = S 1 + £S 2 + ... + £ r ~ 2 S r - 


( 2 . 6 ) 
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where r is any prefixed natural number, Hamiltonian takes the following form 


H = H^ r {j,p x ,ex,e) + e r H r (I, <j),p x , ex, e), 


(2.7) 


77s,r — —p -b Px — 1 + sH-l + ... + £ r 1 77 r _i. 

Here functions Si, 77,: are C°° with respect to I,p x ,ex; Si are continuous functions of 

0- 


Proof. 

Make a canonical change of variables ( I,<f,p x ,x ) h- > ( I,4>,p x ,x ) with generating func¬ 
tions (ESI). The variables are transformed as: 


0 = 0 + £ !?> 

Px = Px + e 2 g, 
££ = ££ + £ 2 7p-. 

opx 


( 2 . 8 ) 


Following the standard procedure of the adiabatic perturbation theory, define functions 
5m such that the Hamiltonian in the new variables does not depend on 0 up to terms of 
order 0(e m ) and S m are continuous and periodic with respect to 0. The procedure is as 
follows: substitute formulas for the change of the “action” and longitudinal variables (p x , x ) 
from EH) to Hamiltonian EH and choose the generating function so that the Hamiltonian 
in the new variables has form EH- 

Thus, one can recurrently define functions S m . For example, an equation for Si looks 
as follows: 


as, p x d(x)d' 

One can choose functions S m in such a way that they have zero average with respect 
to 0. Finally, the Hamiltonian takes the form we need. Functions Si, Hi obviously have all 
the properties mentioned in the theorem. 

Consider the solutions of Hamiltonian equations with Hamiltonians EH and 
7 ~tz,r(J,P,£x,e) corresponding to initial conditions 7(0) = J(0), 0(0) = 0(0), p x (0) = 
p( 0), 7(0) = x(0), where 0 is the phase variable conjugated to variable J. 
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Corollary 2.1 Value of variable I is preserved with an accuracy 0{e r ) on a time interval 
£~ k for any prefixed natural number k. The projection of the ray’s trajectory on p x , ex 
plane on such an interval lies in the 0{e r ) -neighbourhood of curve 7 i.y, )r (J,p,£x,£) = 0. 
Behaviour of variables £x,p x ,(f> is described by the solution of the system with Hamiltonian 
with accuracies 0(£ r ~ k+1 ), 0(£ r ~ k+1 ) , and 0 (£ r ~ k ) respectively (under some natural 
additional conditions given below). 

Proof. 

Similarly to the previous section, consider an auxiliary change of variables with generating 
function 


W = 1(f) + p x X + £S 1 (I,(j),p x ,£x) + . . . + £ r+k 1 S r+k ^i(I,(p,p x ,£x). 

According to Theorem [21 Hamiltonian (12.411 takes the form: 

H = Hz,r+k(i,Px,EX,e) + £ r+k H r+k (I,(j),p x ,£X,£), 

where 0 and x are considered as functions of new variables (I,<f>,p x ,x) . Formulas for the 
change of variables give I—I = 0{e t ), 0—0 = 0(£ r ), x—x = 0(£ r ) and p x —p x = 0(e r+1 ). 
Therefore, it is sufficient to consider behaviour of variables I,(f>,p x ,x. 

Consider an approximate conservation of variable I. The ray’s Hamiltonian in the new 
variables is 


H = Hz,r +k(I,Px,ex,e) + E r+k H r+k {I,<!>,p x ,EX,£). (2.9) 

On each time interval between the reflections of the ray at the upper wall variable / vary 
according to Hamiltonian equation 


j _ ^r+k 9H r+k 

d(j) 

After summing up over all such intervals, it gives an accuracy of conservation 0 (£ r ). 

It follows from J23 that values of variables (p x , x) jump at a reflection of the ray at 
the upper wall. However, the generating function for the change of variables is continuous 
with respect to 0. Therefore, formulas for the change of variables give A p x = 0(e 2 AI) , 
A ex = 0(£ 2 AI ), where A/, A p x , A ex are jumps of variables / , p x . ex after reflection 
of the ray at the upper wall. Using Lemma [2J one estimates value of the jump of variable 
7 as 0 (£ r+k ). The ray reflects at the upper wall not more than 0(s ~ k ) times on the time 
interval £~ k . Therefore, total variation of variable I is 0(£ r ). Thus, the part of the theorem 
concerning variable I is proved. 
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Projection of the rays trajectory on p x , ex plane is in 0(e r ) -neighbourhood of the curve 
'Ht, ,r(J,P,sx,e) = 0 because along the trajectory Hy: t r{Jjpx,£x,£) = 0(e r ). 

Values of variables p,x,ip change continuously, yet values of variables p x ,x, 0 jump at 
a reflection of the ray at the upper wall. It was noted above that values of jumps of variables 
p x ,x,(f) are determined by the value of jump of / at the reflection of the ray at the upper 
wall as A p = 0(e 2 AI) = 0(e r+k+2 ) Ax = O(eAl) = 0(e r+k+1 ) and A0 = O(eAl) = 
0(e r+k+1 ) . Thus, total variations of values of variables (p x ,x,4>) due to reflections are 
0(e r+ 2 ) for variable p x and 0(e r+1 ) for variables x h 0. Total variations of variables 
p x , x , 0 due to reflections are the same. 

Consider the zero approximation of Hamiltonian 7dv,r : 


2 T 2 


0 W F _ 71 I 

s ’ r d 2 {ex ; 


The Hamiltonian system with Hamiltonian r describes the motion of a particle in 


+ p~ 


potential U = ^FTex) ■ Assume (see FigEj) that for the ray under consideration value of 
function F satisfies one of the following conditions: 1) it is larger than value of any local 
maximum of 2 U or 2) the ray can reflect at a potential hump only once or 3) the ray moves 
between two potential humps (resonator). 

Consider case 1). The motion of the ray has a given direction, variables x and x vary 
monotonically. One can consider variables p x and p as functions of variables ex and ex 
respectively. At x — x values of variables p x and p differ by 0(e r ). It takes different time 
for variables x and x to reach a given value x *. This difference can be estimated as follows: 


/ (I -l) dx = 0(0. 

*( 0 ) 

Hence, for a given moment of time t one obtains: ex{t) — ex(t) = 0{e r ~ k+1 ) 1 p x (t) — p{t ) = 
0(e r ~ k+l ). 

The deviation of phases 0 — ip at x = x = x* is : 

AAA _ dAiffifdu + 0(0 = 0(0 (2.io) 

X X J 

where we introduced u = d ^' r . Therefore, for a given moment of time t one gets: 0(t) — 
ip(t) = 0(e r ~ k ). 

In case 2) the ray can change its direction of propagation only once. Here one can consider 
two different regions and use either variable x or p as monotonically changing variable 
in the corresponding region. In case 3) the waveguide is configured as a resonator. One 
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can introduce an “angle” variable in the phase portrait of the system with Hamiltonian 
'H'e, r(I,P,x) as a monotonically changing variable. All arguments here are the same as in 
case 1). 

Corollary 2.2 Value of variable I is preserved with an accuracy 0(e) on a time interval 
e~ r for any prefixed natural number r . 

Corollary 2.3 Behaviour of variables I,(j),p x ,x can be found from the change of variables 
/, (j),p x , x i—► /, (j),p x , x together with approximate expressions defining motion in the system 
with Hamiltonian H^ )r (I,p x ,ex) . An accuracy of such a description is the same as the 
accuracy for variables I,<j),p x ,x according to corolla, ru \2.1\ . 

Remark 4 The Hamiltonian system with Hamiltonian 

v1-2/2 

h =^T+P 2 x -1 ( 2 . 11 ) 

considered on the energy level H = 0 describes the ray’s trajectory with an accuracy 0(e) 
on time intervals of order e _1 . 

Remark 5 Suppose that there is nonhomogeneous medium between the walls such that a 
refraction index n(ex, y) G C°° and trace of the ray in adiabatic approximation is inclined 
to the walls at reflection at a nonzero angle; then Lemma® Theorem® and its corollaries 
10-10 are still valid. 

3 Dynamics of a massive piston in a gas of light particles 

The problem of an adiabatic piston is an important model in statistical mechanics. It is 
considered in the context of attempts to derive thermodynamics laws from the laws of 
mechanics (see, for example, ra.rai- The corresponding system consists of a container with 
a massive cylindrical piston and a gas of identical light particles that move independently 
elastically colliding with the walls of the container and with the piston. 

Let the mass of a particle be equal to 1. Length of the container subtracting thickness 
of the piston L and the number of gas particles are of the order 1. Mass of the piston M is 
large in comparison with mass of the gas. Suppose the piston is at rest at the initial moment 
of time; then energy of the system does not depend on mass of the piston. Hence, one can 
estimate energy of the piston as 0(1) and its typical velocity as 0(-4 =). Therefore, it is 
useful to introduce a small parameter e = -?= . Let us discuss the problem on a time interval 
of order e _1 . 
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Below, to shorten the calculations, we consider the case when there is only one particle 
on each side of the piston (Fig 03). One can easily generalize following arguments for the case 
of any prefixed number of particles. At the end of the section we give the result for the 
general case. Without loss of generality, assume that the particles’ velocities are parallel to 
the axis of the container; thus, the motion is one-dimentional. Let indexes l, r correspond to 
variables of left and right particles respectively. The variables without indexes correspond 
to the piston. Thus, full energy of the system is of the form: 


E = e 


,P 2 


pf 

+ + 


pl_ 
2 ’ 


(3.1) 


where p^ r are momenta of particles, P is momentum of the piston. Denote as xi :T distances 
between particles and the left wall of the container, as X distance between the piston and 
the left wall of the container. The distance for the right particle is calculated subtracting 
the thickness of the piston. 

Suppose the piston is fixed. For each particle make a canonical change of variables from 
(p,x) to “action-angle” variables (7,0) the same way as in Fermi-Ulam problem (Sec. 1). 
Let phases 4>i,r be zero at the walls of the container; then 


vrf, 0G(O, tt), 
vr(2-f), 0G(vr, 2vr), 


0r 


0 e(o, vr), 

*■(2 0€(tt, 2tt). 


(3.2) 


Conjugated “action” variables are /; = , I r = Q . If phases are defined by 

jia, generating functions for the change of variables (p, x ) to “action-angle” have the 
following forms: 


Si{I h x h X) 


S r (I r , x r , X ) 


TTlf, e (0, ?r), 

7rA(2-f), <f>i e (vr, 2 tt), 

n Ir L^x r) ’ & G (°> 7r )> 

7rj r (2 - jjzx), 4> r G {Xi 2vr). 


Make the change of variables with generating function W = PX + Si(Ii,xi,X) + 
S r (I r , x r , X) in the exact system, when the piston is moving. The Hamiltonian in the new 
variables has the form: 


n = e 2 - P 


Y/(0O + L ^ y /( 0r 


7 PI 2 

+ TTTTTr + 


2 t2 


7 PI, 


2X 2 2(L — X) 2 ’ 


(3.3) 


where function /(0) defined by (II .41) . 
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Let us find jump of the “action” variable at an impact. It is natural to fix values of 
all variables except the variables of the colliding particle, obtain the Hamiltonian for the 
particle and consider the conservation law for the value of the Hamiltonian. Let us formulate 
the corresponding lemma for the left particle (the lemma for the right particle is formulated 
analogously). Let the prime denote value of the variable after an impact. 


Lemma 3 This approach defines actual value of the “action’ variable after an impact: 




2PX 2 PX 
7 tM + M + 1 IT M 


(3.4) 


Remark 6 Velocities v and V of the particles of mass 1 and M transform after an 
elastic collision as follows: 


I V' = V + J fc{v-V), 

\ v > = 2V-v + j2 ri (v-V). 

“Action” of the particle and momentum of the piston are Ii 
which follows \3.f\)- 


^ , P = MV, from 


Remark 7 If two particles collide with the piston at the same moment , the further dynam¬ 
ics is not defined. Measure of initial conditions corresponding to such collisions equals 0 . 
Therefore, we do not consider these initial condition. 


Proof. 

An impact is the the passage of (fi through the value tt . Equating values of the energy 
prior to- and after an impact, we find 


7T 2 /f M + 1 _ 7 Tl t ( p I r f (0r ) \ = TT 2 // 2 M + 1 p 4/(0r) 

2X 2 M MX 1 L-X’ 2X 2 M MX 1 L-X 


(3.5) 


After calculations one obtains value of the “action” variable: 


I'i = Ii 


2 X i p + I rf( < Pr 


M + 1 7T 


L-X 


(3.6) 


Let P denote value of momentum of the piston prior to an impact; then 


J'-J 2 X (P l 7Tl1 ) 


(3.7) 
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Finally 




and the lemma is proved. 


2PX 2 PX 

TiM + M + 1 SrM ~~ 


Theorem 3 Values of variables Ii. r are preserved with an accuracy 0(e) on a time interval 


Proof. 

Introduce the normalized momentum of the piston by expression P = eP and consider 
Hamiltonian (13.31) with an accuracy 0(e 2 ): 


H = 


P 


n 2 I? 


2 ' 2X 2 + 2(L-X) 2 £ PL- 


X 


f(4> r ) + 0(e 2 ). 


(3.8) 


Make a canonical change of variables ( Ii. r , <f>i, r , £ _1 P , A") i—> ( I l r . (j)^ , e _1 P, A"), 
eliminating dependence of the Hamiltonian on the phase with an accuracy e 2 . The generating 
function for this change of variables has the form: 

W = -PX + + I r( f r + eStilt, fa P, X) + eS r (I r , <f> r , P, A). (3.9) 


Consider the expressions connecting old and new variables: 


( 






T _ T x ®Sl,r 

71, r — 71, r T e Q J, 1 

7 _ : 8Sl,r 

Vl,r 0i,r A ^ dll r ’ 


P 

X 


p c-2 9S(_ , 2 dS r 

1 " r £ ex T c ax ’ 

y I r 2dSi , 2dSr_ 

A + £ dP +£ w- 


(3.10) 


Substitute these expressions in Hamiltonian (13.81) and choose functions Si >r so that the 
Hamiltonian in the new variables does not depend on the phases in the first order in e. 
Therefore, functions Si, r satisfy equations 


dtp = PXfifa) dSr = P(L- A)/(0 r ) 

d<f>i 7r 2 ’ d(p r 7r 2 


It follows from EH) that functions Si, r are defined up to an arbitrary function of 
variables Ii, r ,P,X. Let this function be equal to 0 and call variables Ii, r the “improved 
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actions”. It is easy to see from (HDD, JHU), that the new phase coincides with the old one 
(( j>i >r = (j)i tr ), and the Hamiltonian has a discontinuity at (fi^ = n . 

Lemma [J| and formulas (EH. El allow to find variation of the “improved action” 
after an impact as 0(e 2 ) . Variation of the “improved actions” between collisions is also 
0(e 2 ) . One can estimate the number of collisions on a time interval e^ 1 as 0(e _1 ). Hence, 
on a time interval e^ 1 “improved actions” vary by 0(e) . The “action” variables are related 
to “improved actions” by formulas for the change of variables Ii, r = Ii, r + £ ^f L ■ Therefore, 
“action” variables vary by 0(e) on the same time interval and the theorem is proved. 


Corollary 3.1 Let 


2 P 2 

H = e 2 — + 


7T 


2 1? 


+ 


7r 


2 j2 


2X 2 2 (L~xy 


Ii r = const, 


(3.12) 


where Ij. r are initial values of the “action” variables. The Hamiltonian system with Hamil¬ 
tonian firm describe the behaviour of variables eP,X with an accuracy 0(e) on a time 
interval e^ 1 . 


One can similarly consider the system with any prefixed number of particles. Value of 
the “action” variable of each particle is preserved with an accuracy 0(e) on a time interval 
e -1 . The Hamiltonian system with Hamiltonian 


H = e 2 


P 2 


+ 


7T 


-T 2 + 

, 1 T 


7r 


(3.13) 


2 M ' 2X 2 EJ 2(L — X) 2 E ’” 
describes behaviour of variables eP,X with an accuracy 0(e) on the same time interval. 
Here /1 r and /1 1 are the sums of squared initial values of “action” variables for the left 
particles and for the right particles respectively. (This result was first derived in |2] by a 
different method.) 

Remark 8 In the considered approximation the piston oscillates in potential U = + 

2 (L-x)^ kr (Figtity- 
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Phc. 1: A ball between slowly moving walls. Phc. 2: A phase portrait for the motion of the 

particle. 



Phc. 3: A planar slowly irregular waveguide. 
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Phc. 4: The plot of function 2 U and levels of 
function F. 



Phc. 5: The system of a piston and two particles. 


Phc. 6: The effective potential energy of the pis¬ 
ton. 
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